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("^ ' Abstract. - Extending finite size scaling theory to the many body ground state, one finds 

0^ I that Coulomb repulsion can drive a system of spinless fermions in a random potential from the 

^\ , Anderson insulator (Fermi glass of localized states) towards a new extended phase in dimension 

d = 2. The transition occurs at a Coulomb energy to Fermi energy ratio r^ ~ 4, where a change 

C^ ' in the nature of the persistent currents has been previously observed. Relevance to the recently 

C^ I observed 2d metallic phase is suggested. 

o 
o 



Following the scaling theory of localization y, all 2d systems of electrons (or holes) are 
localized when electron-electron interaction is negligible. This absence of metals in two 
dimensions is nowadays questioned after the observation of a transition towards a metallic 
behavior in silicon based 2d electron gases [Q, H, ^ when the carrier density is varied. This 
unexpected phenomenon occurs also for hole gases in GaAs heterostructures |^, ||, SiGe [Q 
and InAs quantum wells S. The transition is observed at very low carrier concentrations, the 
charge spacing being > lO^A in certain cases p]. T his suggests that Coulomb repulsion is the 
driving mechanism for the transition. In ref. [|10|, a very clean heterostructure was studied, 
and the transition was observed at Tj « 35, close to r^ « 37 for which Wigner crystallization is 
expected without disorder. In more disordered samples, the transition occurs typically around 
Tg Ri 10 and kpl w 1, kp denoting the Fermi wave vector and I the elastic mean free path. It is 
also at Ts « 10 that charge crystallization is numerically observed in small disordered clusters 
[pT, 12 1 . This gives arguments to associate the observed transition to the quantum melting of 



a kind of pinned Wigner crystal. In this case, the metallic phase should cease to exist at a 
weaker rg also and the re-entry towards a Fermi glass of weakly interacting localized particles 
should occur. The re-entry has been observed in Ref. at r^ w 6 (see also Ref. H). Those 
regimes are out of reach of the Landau theory of disordered metals O, H valid for large kpl 
and small rg. 

The interplay between disorder and electron-electron interactions is attracting increasing 
theoretical interest [n5|, though many studies are devoted to large energy excitations |16| of 
a few particle states instead of the ground state properties. In Ref J12J, it was shown that 
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the ground state of a system of spinless fermions goes from an insulating phase (Fermi glass, 
Ts < rf) towards another insulating phase (pinned Wigner crystal, rg > r^) through a 
new intermediate quantum phase. This conclusion was drawn from a study of the persistent 
currents supported by the ground state in small clusters. A 2d torus geometry enclosing a flux 
(j) around the longitudinal direction was considered. The ratio rf characterizes the suppression 
of the transverse current It (not enclosing 0) . Below rf , the pattern of driven currents has a 
2d topology and the longitudinal current // (enclosing 0) can be paramagnetic or diamagnetic, 
depending on the sample. Above rf , the flow pattern has an ordered Id topology [|l2|, |lj] and 
the sign of /; becomes sample independent. The higher ratio r^ characterizes the suppression 
of/; and charge crystallization in a random substrate. For 0.3 < kpl < 3, one obtains r^ ~ 10 
and rf « 4, in agreement with ratios r^ « 10 where the metal-insulator transition is observed, 
and with a ratio r^ « 6 where the re-entry has been reported in Ref. P| . Though small clusters 
exhibit an enhancement of /; for intermediate r^ , exact diagonalization does not allow to vary 
system size and to establish that the intermediate phase is metallic at the thermodynamic 
limit. To understand how the cooperative Coulomb behavior begins to destroy Anderson 
localization of weakly interacting particles, and if it drives the system towards a delocalized 
phase befo re y ielding charge crystallization, we need to extend the finite size scaling method 
[ |l8[ |l9| , po| , |2l[] successful for describing single particle delocalization to the many body ground 
state. We do not study in this work large Coulomb repulsions where charge crystallization sets 
in, but using approximations valid for weaker repulsions, one finds that the Fermi glass melts 
at Ts « rf" ~ 4 (in good agreement with Ref. |1^) to give rise to a more fiuid phase for the 
particles. Using a low carrier density, we characterize this melting by a suitable localization 
length which displays in d = 2 a scaling behavior similar to the one body (IB) localization 
length in d = 3 at the mobility edge, i.e. more generally a behavior characteristic of second 
order phase transitions. This confirms that spinless fermions have an intermediate delocalized 
phase between the Fermi glass and the pinned Wigner crystal. 

Anderson IB localization is a complex phenomenon which is analytically tractable in a few 
limits: quasi-one dimension [g2|, |2^ and Bethe lattice [Q. For dimensions d = 2 and 3, our 
knowledge of IB localization is mainly based on numerical works [|l8[ |l^, ^, ^ . The successful 
method consists in evaluating the IB localization length Li{L) of a system of finite size L and 
to verify a scaling ansatz inspired |Iq] from the theory of second order phase transitions: 

The lengths Li (L) calculated for different system parameters can be mapped onto a universal 
curve Fd assuming a single scaling length Li{d) which characterizes the d-dimensional lattice. 
For the IB problem, it was convenient to consider d-dimensional strips of finite transverse 
section L'^~^ . Li{L) was defined as the inverse of the smallest positive Lyapunov exponent of 
the appropriate product of transfer matrices, which is a self averaging quantity. The ansatz (|lh 
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was verified in dimensions d> I, and a metal- insulator transition was obtained [|T8[ |l^, ^, 
for d = 3 with Li{3d) diverging at the mobility edge. In two dimensions, large L studies |20 
led to the conclusion that Li(2d) diverges only in the clean limit, in agreement with Ref. |]l[ . 

We extend this finite size scaling method to the ground state of N spinless fermions with 
Coulomb repulsion in a random potential defined on a square lattice with L^ sites. The 
Hamiltonian reads: 

H = -tY^clc,+Y^v.n, + uY^'^. (2) 

c\ [ci) creates (destroys) a particle in the site i = {ix,iy), i is the strength of the hopping terms 
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Fig. 1. - A: IB density of states pi{E) and number Mi k, L\ of occupied sites (insert) as a 
function of the energy for W = IQ and L = 31. B: Convergence of the number Nsd of SD 
occupied by the ground state when Nh increases, i = 31 (upper curve) and L — 28 (lower 
curve) at t/ = 1.5. C: Variation of Nsd as a function of U for L ~ 24 (circle), 28 (square) and 
31 (diamond). D: Number Mg of occupied sites per particle as a function of U for the three 
values of L (same symbols than beside). 



between nearest neighbors (kinetic energy) and r^ is the inter-particle distance for a 2d torus 
(periodic boundary conditions) . The random potential Vi of the site i with occupation number 
Tii = c\ci is taken from a box distribution of width W . The interaction strength U yields a 
Coulomb energy to Fermi energy ratio r^ = U / {2t^/Tm2) for a filling factor rig = N jl? . 

To investigate the melting of the Fermi glass by Coulomb repulsions, one begins with a 
low density of occupied IB states such that the charge spacing is larger than L\ without 
interaction. The size L varies from L = 24 (iV = 3) to L = 31 (N = 5), TV being chosen to 
keep rig — 1/192. An intermediate size L — 28 (N — 4) with almost the same density (1/196) 
is also considered. 

Let us first discuss some features of the IB spectrum. If we take a small value for W, strong 
lattice effects remain for L w 31, complicating the extrapolation towards the thermodynamic 
limit. The lattice effects disappear at larger W, but another difficulty remains, due to the low 
density rig numerically accessible. The low energy tail of the IB spectrum is made of impurity 



4 EUROPHYSICS LETTERS 

states where the particles are simply trapped at some site i of exceptionally low Vi. Their 
energies are below the band [— 4i, 4i] of the clean system. This is a trivial localization with 
Li < 1 and irregular energy spacings given by the tails of the distribution of the Vi. We are 
not interested to study the detrapping of those impurity states by Coulomb repulsions, but by 
the delocalization of genuine Anderson localized states, of energy E > — 4i. The localization 
of those states results from the multiple scattering of plane waves, while the impurity states 
involve mainly a few evanescent waves. With a higher density, it would be easy to fill the tail 
of IB band and to put the Fermi energy inside the band of Anderson levels of energy > — 4t, 
having regular energy spacings Ai oc 1/L^ and similar Li (see Fig. 1 A). With a low density, 
in order to observe nevertheless a delocalization resulting from the mixing by the interaction 
of SD built out from Anderson localized states (say with Li ~ 4), we calculate the L^ IB levels 
for W = 10 and we ignore the L^/2 first IB levels, considering only as possible IB states the 
L'^/2 remaining ones. Doing so, the states near the Fermi energy of the non interacting system 
are not simple impurity states. We have checked that all the results presented hereafter are 
identical if we get rid of a smaller fraction of IB states, as far as the Fermi energy at [/ = is 
located in the constant part of the IB density of states (Fig.l A). For W^ = 10 and IB states 
around the IB band center, one can roughly estimate that kpl « 1. From this restricted subset 
of IB Anderson localized states, the Nh first SD ordered by increasing energy are calculated. 
We assume that such a procedure captures genuine 2d ground state physics, since the low 
energy excitations begin with a (large) IB energy spacing Ai, and not with a (very small) NB 
energy spacing Ajv oc L~^^ . 

The sizes Mh = (i^)!/(iV!(L^ — N)\) of the corresponding Hilbert space are huge and only 
the projection of the Hamiltonian (0) onto smaller subspaces can be diagonalized (Lanczos 
method). For weak [/, we consider the subspace spanned by the Nh Slater determinants 
(SD) corresponding to the first low energy states of the non interacting problem. The N 
body (NB) ground state |^o > is obtained after diagonalizing the truncated Hamiltonian. 
1*0 >= Ej C^\SD^ > is typically extended in the SD basis over Nsd =< Ej \Ci\'^ >^^ SD, 
the brackets denoting the average over 5 x 10^ disordered samples. In Figs. 1 B-C, one can see 
that Nsd remains stable when Nh > 500 at a value negligible compared to Nh = 1000 (value 
for Nh assumed hereafter) . Though the stability of our results have been checked when Nh 
varies inside a range negligible compared to Mh , making difficult to rule out the existence of 
a slow variation which cannot be detected inside a too narrow range, we assume that |^o > 
has a negligible chance to have a significant projection outside the 1000 first SD when i < 31 
and U < 1.5. 

Figure 1 D gives the number AIs — N < J^iPf >^^ o^ occupied sites per particle in the 
ground state, pi —< *o|^i|*o > denoting the charge density of the ground state at site i. 
For U = 0, Ms is obviously smaller than Lf, a limit where the N occupied IB states do not 
overlap. When one goes towards the thermodynamic limit, Ms is almost size independent for 
the Fermi glass {U < 1), but varies as a function of L when U > 1. This provides a first 
evidence that one has a delocalization effect at an interaction threshold C/ « 1. The glassy 
ground state melts towards a more fiuid phase which fills a larger fraction of the sample. This 
observed delocalization cannot be related to charge crystallization, since Ms increases while 
Ms -^ 1 when U — > oo. 

To characterize the NB ground state by an appropriate localization length, we consider 
the change Spj of the charge density induced by a small change Svi of the random potential 
Vi located at a distance r = \i — j\. For a Slater determinant made with N occupied IB 
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Fig. 2. - A: < ln|(5p(r)| > as a function of r for L = 31 at [/ = (circle), 0.25 (square), 0.5 
(diamond), 0.75 (triangle up), 1 (triangle left), 1.25 (triangle down) and 1.5 (triangle right). 
The lengths £,l{U) are given by the slopes of the thick lines. B: (_l as a function of U for 
L = 2A (circle), 28 (square) and 31 (diamond). 



eigenfunctions ( ■(/'q)j fi^'st order perturbation theory gives: 



N 

5p.j = 2Sv^ Y^ Y^ 

a=l f3=ia 



E[3 — Ea 



2r 

oc exp , 



(3) 



the index /3 varying over the IB spectrum. Therefore, the change Sp remains localized 
on a scale ^l ~ Li/2 without interaction. We study how ^^ depends on U for different 
values of L. To improve the statistical convergence, we calculate more precisely the change 
Sp(r) = ^ ^Pr,jy of the charge density on the L sites of coordinate jx = r yielded by 
the change wo,i ~^ "^0,2 (1 + 1/100) for the L random potentials of coordinate i^ = 0. The 
quantity Sp{r) for M^ = 10 has a broad symmetric distribution of amplitude |(5p(r)| which is 
almost log-normal for U = 0. There is no underlying law of large numbers which tells us what 
is the right self-averaging scaling variable, as Oseledec's theorem for the quasi-ld IB problem 
[ p^Sl . Though the log-normal character of the distribution becomes less pronounced when U 
increases, it still makes sense to characterize the typical strength of the fluctuations by 



T 

\Sp{'^)\typ = exp < ln\Sp{r)\ >cx exp--—. 



(4) 



where the brackets denote the average over 5 x 10^ disordered samples. We extract the range 
^L of the perturbation from the r dependence of this typical value inside a square of size L. 

An exponential decay of |(5p(?')|typ occurs only over a scale r << L/2 since the boundary 
conditions are periodic. In Fig. 2 A, one can see how the length ^/, are obtained from the 
slope of the linear parts of the curves (thick straight lines). Fig. 2 B gives how ^^ depends on 
U for the three considered sizes. This figure conveys a very different information than Fig. 1 
C. The number Nsd of SD over which the NB ground state has important projections linearly 
increases as a function of U, without exhibiting a change of behavior. This rules out that U 
can induce a sharp localized-delocalized transition in Hilbert space of the type suggested in 
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Fig. 3. - A: Localization lengths S,l{U) as a function of L for U <1 (filled symbols): (triangle 
up), 0.25 (triangle left), 0.5 (triangle down), 0.75 (triangle right), 0.85 (circle) and 1 (square) 
and U > \ (empty symbols): 1.15 (square), 1.25 (triangle), 1.35 (diamond) and 1.5 (circle). 
B: Ratios ^l{U)/L mapped onto the scaHng curve F2 as a function of ^{2d)/L. The two 
dimensional scaling length ^{2d) is given in the insert. 



Ref. pq | for metallic quantum dots. The behavior of £,l{U) is much more interesting. When 
L = 31 and U > I, the range of a local perturbation has a sharper increase which can only be 
due to the nature of the last SD participating to the NB ground state, and not to their number 
NsD- The last SD participating to the ground state when U > 1 should not be made of IB 
states localized in the vicinity of those playing a role when U < 1. This may indicate that the 
underlying delocalization mechanism is related to variable range hopping between localized IB 
states yielded by Coulomb repulsions. This possibility was discussed in Rcfs. [^6|, ^, psj . 

The scaling analysis is shown in Fig. 3. The size dependence of £,l is presented in Fig. 3 A 
for increasing Coulomb repulsions. One finds the behavior typical of a transition: ^l converges 
towards a finite value when U < Uc (localized Fermi glass); diverges linearly as a function of L 
at U = Uc ^ 1 (critical point) and diverges faster than linearly when U > Uc (extended phase). 
This is exactly the behavior |18) which characterizes the Anderson transition at d = 3 for the 
IB spectrum. In Fig. 3 B, one verifies the scaling ansatz (|l|), where £,l and (,{2d) play the 
role of Li{L) and Li{d) respectively. All the data of Fig. 3 A can be mapped onto a universal 
curve F2 shown in Fig. 3 B, assuming the scaling length ^{2d) given in the insert. When 
U < Uc, this length characterizes the localization of the effect of a local perturbation of the 
substrate in the two dimensional thermodynamic limit. Considering that the delocalization 
threshold is located at Uc ^ I, a. power fit of (,{2d) (x \U — Uc\~^ yields a rough estimate 
for the critical exponent v k, A. More detailed and accurate studies of the vicinity of the 
delocalization threshold are necessary for confirming this value for v. Very often, additional 
corrections ex L~°' to the scaling ansatz occurs for small sizes. We point out that our results 
can be fitted by a simple linear law ^l — 0.17i for U ~ Uc- This is a further indication 
that the simple ansatz (M describes scaling for L > 24, without noticeable additional L~° 
corrections. The obtained interaction threshold Uc ^ 1 gives r^ w 4, which agrees with the rf 
given in Ref. |12 for similar values of kpl and is consistent with the threshold r^ « 6 where 
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the re-entry is reported in Ref . g . 

In summary, we have shown that there is a mching of the Fermi glass towards a new 
delocahzed phase at Tg = rf « 4 and kpl ss 1. In the glassy state, the particles occupy a tiny 
fraction < TigL^ of the sample and the effect of a local perturbation remains localized. Above 
rf , the ground state becomes extended over the whole sample, as shown by the divergence of 
the range £,{2d) characterizing the two dimensional thermodynamic limit when r^ —> rf . The 
range ^x, calculated in a square of size L satisfies a finite size scaling ansatz consistent with a 
second order quantum phase transition and a power law divergence of £,{2d) at rf . The effect 
of a local perturbation (motion of a single atom) have interesting implications for 1// noise 
in the metallic p3| and localized [pffl regimes. Therefore, the implications of our numerical 
results might be checked by suitable noise measurements in low density 2d gases of electrons 
or holes. 

This work is partially supported by a TMR network of the EU and by a Franco-Israeli 
AFIRST grant. 
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